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We study 2D wedge wetting using a continuum interfacial 

Hamiltonian model which is solved by transfer-matrix 
methods. For arbitrary binding potentials, we are able to 
exactly calculate the wedge free-energy and interface height 
distribution function and, thus, can completely classify all 
types of critical behaviour. We show that critical filling is 
characterized by strongly universal fluctuation dominated 
critical exponents, whilst complete filling is determined by 

the geometry rather than fluctuation effects. Related 
phenomena for interface depinning from defect lines in the 
bulk are also considered. 

PACS numbers: 68.45. Gd, 68.35. Rh, 68.45.-v 



At present, experimental methods allow the shape of 
solid surfaces to be controlled at a nanoscopic level jij. 
Fluids confined by such structured substrates can 
exhibit quite distinct adsorption characteristics 
compared to that occuring for planar systems || . This 
includes new types of interfacial phase transitions and 
critical phenomena which are not only of fundamental 

interest but may well play an important role in 
developing technologies such as super-repellent surfaces 
H , self assembly of three-dimensional structures 0] or 
micro-fluidics 0, among others ||. An interesting 
example of these phenomena which is recently 
attracting new interest is the so called filling or wedge 
wetting transition of a fluid adsorbed in a wedge |7fl§] , 
formed by the junction of two flat walls tilted at angles 
±a to the horizontal as shown in Fig. 1. 
Thermodynamic arguments predict that a wedge-gas 
interface is completely filled by a liquid phase (at bulk 
liquid-gas coexistence) for temperatures T > T Q , where 
the filling temperature T a is lower than the wetting 
temperature T w of the planar (a = 0) wall In fact, 
according to these macroscopic arguments, the location 
of the filling transition phase boundary is beautifully 
expressed in terms of the contact angle O^CH of the 
liquid drop on the planar substrate 



e w (T) 



(i) 



Thus, the liquid completely wets the wedge when the 
contact angle is smaller than the tilted angle a. 
Interestingly, this macroscopic result was predicted and 
confirmed experimentally |l0| eight years before the 
seminal paper by Cahn on wetting in planar surfaces 

in. 



Recently, the macroscopic prediction (Q) has been 
supported by mean-field analysis of a model system 
which also suggest that the qualitative order of the 
filling transition (first-order or continuous) follows that 
of the planar wetting transition Q). Thus, for planar 
substrates exhibiting critical wetting transitions, the 
wedge offers two new examples of interfacial-like critical 
phenomena in which the interface height £q (measured 
from the bottom of the wedge) diverges as the 
temperature and chemical potential are varied. 
Borrowing from the vocabulary used for wetting, we 
refer to the wedge filling transition ocurring as T — > T~ 
(at bulk coexistence) as critical filling. In contrast, by 
complete filling, we refer to the divergence of £q for 
temperatures T > T a as the bulk chemical potential /i is 
increased towards saturation ^ sa t(T). However, apart 
from a few limited results available for the corner 
wetting transitions (corresponding to a = 7r/4 and 
restricted to short range forces) 12 1, there has been no 
discussion in the literature of fluctuation effects, scaling 
regimes and universality classes for such filling 
transitions and how these compare with the rich 
phenomenology known for wetting ]l3| . With this aim 

in mind, we have studied an effective interfacial 
Hamiltonian model of filling in (bulk) dimension d = 2 
and derived exact elegant results for various quantities 
of interest such as the excess wedge free-energy and the 
probability distribution for the interface height. The 
formal analysis can be carried through for arbitrary 
choices of binding potential (i.e., all range of forces) 
allowing a complete classification of the critical 
behaviour and the identification of universality classes. 
We will show that both critical and complete filling 

transitions are characterised by universal critical 
exponents independent of the intermolecular forces 
(unless they are unphysically long-ranged). 
Interestingly, whilst the universality of the critical filling 
transition arises directly as a consequence of strong 
fluctuation effects, the universality encountered at 
complete filling has a geometrical and thermodynamic 
origin which is correctly captured by mean-field theory. 
In addition to this, we are able to verify the validity of 
the thermodynamic prediction (Q) and also predict very 
similar phenomena occurring when interfaces are pinned 

in the bulk along a bent line of weakened bonds. 
To begin, we introduce our model. We denote the local 

height of the liquid-gas interface relative to the 
horizontal by y(x) and the height of the wall itself by 
z(x) = a\x\. The relative local height between the two 
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is written £{x) = y{x) — z(x). The horizontal dimension 
of the wall/interface extends over the range 
[-X/2, X/2] and periodic boundary conditions are 
applied at the ends; i.e., all configurations satisfy 
y(X/2) = y{-X/2) (see Fig. (1)). The interface 
interacts with the wall via a binding potential arising 
due to intermolecular wall-fluid and fluid-fluid forces 
and is also subject to thermal fluctuations governed by 

the stiffness E. We will concentrate on continuum 
wall-fluid systems and identify the stiffness E with the 
surface tension. We also set ksT= \ for convenience. 
For rather open edges, it is already known from earlier 
mean-field studies that the filling transitions are well 
described by the simple interfacial model 



H[y] 
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which can be justified from analysis of a more general 
non- linear model in the small angle a limit Q . Before 
we outline our calculation and present our main results, 

we make some preliminary remarks which serve to 
establish our notation and also provide some points of 
comparison. These concern the planar limit a — 0. For 
this case, it is well known that the partition function 
Z^(£i, £2; X) corresponding to the planar fluctuation 
sum with fixed boundary conditions £(—X/2)=£i and 
£{X/2)=t2 is given by the spectral sum (or integral if 
scattering states are present) Q 



Z V (£ 1 ,£ 2 ;X) = J2 e 



E n X 



(3) 



n=0 



where the eigenfunctions and eigenvalues satisfy the 
Schrodinger equation 



f 
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Thus, in the thermodynamic limit X — * 00, the excess 
free-energy is given by Eq which is in turn related to 

the contact angle (defined only for T < T w and 
fx = Hsat(T)) by E — — E0J/2, valid for small angles 
Qtt described by the interfacial model. Similarly, the 
normalised probability distribution for the interface 
height is P v {£) = \^)q{1)\ 2 . For later purposes, it is also 
convenient to define the matrix elements 



m\f(£)\n)= d£ Vm(*)/Wn(*)- 



(5) 



These will appear in our solutions for the wedge 
free-energy and probability distribution function. The 
binding potentials that we consider are of the usual 
form (TJ 



W{£) = hl+ + 



£>0, 



(6) 



with h cx (p sa t -/j), fl and b are effective Hamaker 
constants and q>p>0 allow for general types of 



intermolecular potentials. Now, recall that for two 
dimensional critical wetting transitions, the critical 
behaviour gencrically belongs to one of three scaling 
regimes depending on the values of p and q. Specifically, 
for q>p>2, p<2 but q>2 and, finally, p<q<2, the 
behaviour falls into the strong, weak and mean-field 
fluctuation regimes with true universality only 
characteristic of the former [jl3| . For example, the mean 

interface height (£} v ~ (T w — T) _/3s diverges with 
critical exponents f3 s = 1, f3 s = 1 j{2 — p), f3 s — 1 /{q — p) 
in the three regimes, respectively. For complete wetting, 
corresponding to h — * for T > T w , there are only the 
weak and mean-field fluctuation regimes and the 
interface height diverges as {£} 7T ~ hr 1 / 3, and 
(£) 7T ~ h~ x /( p+1 \ respectively. These remarks will serve 

to illustrate how different the critical behaviour at 
wedge filling is to standard planar wetting transitions. 
We now turn our attention to the wedge geometry and 
outline the calculation of the partition function 
Z^ edge (X) for the present periodic system. The 
advantage of this choice of boundary conditions is that 
it allows us to extract the excess wedge free energy 
F W edge(ct) rather easily from 
FLd 9 e(X) = -\nZP edg? (X). To see this, note that in 
the thermodynamic limit X — * 00, the periodic system 
reduces to two independent wedges because one must 
also consider the contribution from the (inverted) wedge 
at x~±X/2 characterised by an angle —a. First, from 
F^ edge (X) we subtract the free-energy F% (X) of a 
planar system (with periodic boundary conditions) 
extending on the same area, i.e., a flat wall tilted an 
angle a to the horizontal. This can be easily calculated 
with the same Hamiltonian, Eq. (^), but with z(x)~ax 
(for all x in the range —X/2<x<X/2). This defines 
the excess periodic wedge free-energy: 



wedge 



(X)-FP(X). 



(7) 



Thus, in the thermodynamic limit, we can write 
lim AF^ ed (X) = F wedge (a) - F wedge (-a), (8) 

X — >oo * 



illustrating the independent contributions from the 
wedge and inverted wedge. The partition function 
Z^ edge (X) is given by the fluctuation sum over all 
graphs y{x) or, equivalently, over all relative positions 
£(x) = y(x) — z(x). Making this change of variable, we 
can rewrite the energy, Eq. (Q), of a configuration as 



H[£] 



-a 2 X 
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dx 

-X/2 
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(9) 



where £ a = y(0) and £ e = y(\X/2\) - z(\X/2\) are the 
midpoint and edge (relative) interface heights. 
Consequently, 
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traced to the asymptotics of the planar ground state 
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Substituting the quantum mechanical result, Eq. (||), 
and taking the thermodynamic limit, we arrive at the 
general formula for the wedge free-energy, valid for all 
binding potentials, 



F wedge (a) = -ln(0\e^ al \0} 



(11) 



where the inner product is defined in terms of the usual 
planar system eigenfunctions (see Eq. (|5|)). Proceeding 

in this way, it is also possible to calculate the 
probability distribution V{£\ x) for finding the interface 
at height £ from the wall at a distance x along it. We 
omit details and simply quote our final result obtained 
in the thermodynamic limit Jlq] , 



|0) M(t)MZ) e-( £ "- B °)N 
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Note that, since the thermodynamic limit is taken first 
(at finite x), this result pertains to a single wedge 
system. There is no contribution due to the inverted 
wedge at x — ±X/2. For the probability distribution of 
the inverted wedge, one simply reverses the sign of a 
and replace \x\ by \x — X/2\. Note that when a = 0, the 
wedge free-energy vanishes and the probability 
distribution reduces to the standard planar result 
Vtt(£). In addition, for a^O, only the n = term 
survives in the limit X — > oo so that V[£\ x) — > V n (£) 
infinitely far from the wedge bottom. The expression 
for the probability distribution simplifies considerably if 
we consider the local height probability at the 
mid-point x = 0. Writing V(l ) = V{£; 0), we find 



V{£ ) 



(0 |e 
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which is one of the central results of our paper. We note 
that both Eq. (|ll|) and Eq. ( p^| ) are consistent since 
from the Hamiltonian definition the mean mid-point 
interface height satisfy 2S(4) = —dF we d ge (a) /da. It is 
also possible to expand V(£;x) for small x and derive 
further explicit results which allow the calculation of 
the curvature of various local operators at x — 0. In the 
present paper, however, we simply concentrate on the 
properties of the wedge free-en ergy F we d ge (a) and 
mid-point height distribution p6| . The essential 
observation to make here is that, relative to the planar 
distribution function Vn(£), the mid-point height 
probability V(£o) has an exponential boost factor 

which decreases the pinning effect of the binding 
potential (provided a > 0, of course). Due to this 
exponential term, the location (phase boundary) and 
character of the filling transition can immediately 
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wave function ipo(£)- If the decay of this function is too 
slow, the wedge distribution V{£q) is no longer defined 
and the wedge is filled with liquid. Now, for bulk 
coexistence (/i = 0) and subwetting temperature 
(T<T W ), the asymptotic decay of ipo(£) has the same 
functional form for all potentials of the form (H): 



,-Ee,(T) t. 



h = 0, T<T W , l^oo 



(14) 



where the specific p and q dependence only enters 
implicitly through the temperature dependence of the 
contact angle T (T) = \j2\E a \/Y>. Thus, the location of 
the wedge filling transition within the present model 
exactly matches with the thermodynamic prediction, 
Eq. (Q). Moreover, since 6 W (T) is analytic away from 
T w , we may identify \Q V (T) — a\ oc T a — T and derive 
the universal critical singularities for critical filling 



F we dge{a) ^ In (T a - T), 



<M ~ (T a -TY 



(15) 



valid for all intermolecular forces (provided p> 1). 
Therefore, the universality of the critical wedge filling 
transition far exceeds that encountered at critical 
wetting, and is ubiquitous to all realistic solid-fluid 

interfaces. We also note that the transition is 
fluctuation dominated since the mid-point roughness 
= yj (£% ) — (£a ) 2 also diverges with the same power 
law. The universality of the critical exponents for 
critical filling is one of the central predictions of our 
paper, and can be traced to the large scale interfacial 
fluctuations occuring at the bottom of the wedge. In 
contrast, mean-field calculations which ignore 
fluctuations effects, predict highly non-universal critical 
behaviour for critical filling. For example, minimization 
of the Hamiltonian, Eq. (||), leads to the prediction that 
the mid-point height diverges like (£q) ~ (T a — T) _1 / p 
Jl6| . In d=2, this prediction is only valid for systems 
with p < 1 which do not correspond to any known 
physical forces. 
Next, we turn our attention to complete filling ocurring 
for h — > and T>T a . For h ^ 0, the asymptotic decay 
of ipo(£) is faster than exponential due to the divergent 
linear term in W{£). This dominates the mid-point 
probability function V{£q) at large distances and 
ensures that the wedge is only partially filled when the 
system is out of bulk two-phase coexistence. The 
singularity arising in the evaluation of the wedge 
free-energy and mean height (^o) can be calculated 
using standard techniques. For the mean height (£o), we 
find that the leading order behaviour as h — > is 

(to) = < (16) 



so: 

2h 



T>T lL 



3 



which is again independent of the intermolecular 
potential exponents p and q. In fact, these results are 
identical to those derived in the mean-field analysis of 
K. Rejmer et. al. equivalent to a simple minimization 
of the Hamiltonian (||). Thus, the universality of the 

complete filling exponents has a geometrical and 
thermodynamic origin rather than being a fluctuation 
related effect. We mention here that the prediction ( [l6| ) 

for the case T > T w is consistent with earlier 
solid-on-solid model calculations of complete wetting at 
a corner . In the light of these transfer-matrix and 

mean-field results, we conjecture that the critical 
singularities occuring at complete filling are of the form 
(16), independent of the dimensionality. 
To end our paper, we mention that very similar 
phenomena also occur for interface pinning in the bulk. 
Recall that an interface is always pinned along a 
straight line of weakened bonds, with dcpinning only 
occuring as the strength of the bonds approaches the 
bulk value |fl3fl . However, this is not the case if the line 
of weakened bonds has a bend in it, as can be seen from 
the present transfer-matrix analysis. To model this 
system, we use the same function z(x) = a\x\ to 
describe the local height of the line of weakened bonds 
but with a square-well potential of depth U and range 
R/2, chosen to mimic the local energy cost due to the 
bond weakness (restricting ourself to systems with 
short-ranged forces). Whilst in the planar system 
(a = 0) the interface has equal probability of being 
found above and below the line, the bend at x = 
breaks the symmetry and significantly enhances the 
probability of finding the interface above the line of 
weakened bonds. In fact, the interface unbinds and 
depins from the defect line at a non-zero value of the 

weakness parameter U a satisfying Eo(U a ) = £ev 2 /2 
where Eq is the ground state energy of the square-well 
potential trivially found from solution of Eq. (^J). 
In summary, we have shown through exact 
transfer-matrix calculation, that critical filling of liquid 
in a 2D wedge is characterized by strongly universal 
critical behaviour. Whilst critical filling is dominated 
by fluctuations, these do not affect complete filling 
which depends solely on the system geometry and 
whose critical behaviour is correctly captured by 
mean-field theory. Similar behaviour is expected for the 
3D wedge. We also point out similar behaviour occuring 
at defect lines in the bulk. 
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FIG. 1. Schematic illustration of an interface configuration 
y(x) in the wedge geometry. The function £(x) denotes the 
local distance to the wall. Other notation is defined in the 
text. 
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